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The Approximate Capacity 
of the MIMO Relay Channel 

Xianglan Jin and Young-Han Kim 


Abstract —Capacity bounds are studied for the multiple- 
antenna complex Gaussian relay channel with ti transmitting 
antennas at the sender, r2 receiving and t2 transmitting antennas 
at the relay, and rs receiving antennas at the receiver. It is shown 
that the partial decode-forward coding scheme achieves within 
min(ti,r2) bits from the cutset bound and at least one half of 
the cutset bound, establishing a good approximate expression of 
the capacity. A similar additive gap of min(ti +t2, ra) + r2 bits is 
shown to be achieved by the compress-forward coding scheme. 


I. Introduction 

The relay channel, whereby point-to-point communication 
between a sender and a receiver is aided by a relay, is a 
canonical building block for cooperative wireless communi¬ 
cation. Introduced by van der Meulen ||T1, this channel model 
has been studied extensively in the literature, including the 
now classical paper by Cover and El Gamal IS). The problem 
of characterizing the capacity in a computable expression, 
however, remains open even for simple channel models, and 
consequently a large body of the literature has been devoted 
to the study of upper and lower bounds on the capacity. 
Reminiscent of the max-flow min-cut theorem El, the cutset 
bound was established by Cover and El Gamal ||2l, which 
sets an intuitive upper limit on the capacity. On the other 
direction, there are a myriad of coding schemes, typically 
referred to as “*-forward” a, each establishing a lower 
bound on the capacity. Among these, the two most versatile 
coding schemes are partial decode-forward a Th. 7] and 
compress-forward El Th. 6], which are complementary to 
each other (providing digital-to-digital and analog-to-digital 
relays, respectively) and have been successfully extended to 
general relay networks for unicast, multicast, broadcast, and 
multiple access ||5|, fb), Q, lH. 

The Gaussian relay channel, whereby the signals from the 
sender and the relay are corrupted by additive white Gaussian 
noise, is one of the most basic channel models studied in 
the literature. The capacity of the Gaussian relay channel, 
however, is again unknown for any nondegenerate channel 
parameter. Instead, the following results have been established 
for single-antenna Gaussian relay channels. 
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• Partial decode-forward and compress-forward, respec¬ 
tively, achieve within one bit from the cutset bound 191 , 

Co). 

• Partial decode-forward, which is a superposition of 
decode-forward and direct transmission, reduces to the 
better of the two M. 

These results establish simple approximate expressions of the 
capacity, which are particularly useful in high signal-to-noise 
ratio (SNR). A natural question arises on how these results can 
be extended to multiple-antenna (also known as multiple-input 
multiple-output or MIMO) Gaussian relay channels. 

Capacity bounds for MIMO relay channels have been 
studied in several papers. By convex optimization techniques 
EH, Wang, Zhang, and Hpst-Madsen ifT^ derived upper and 
lower bounds based on looser versions of the cutset bound and 
the decode-forward bound. These results have been improved 
by more advanced coding schemes (partial decode-forward 
and compress-forward) with suboptimal decoding rules by 
Simoens, Munoz-Medina, Vidal, and del Coso M and Ng 
and Eoschini na. The usual focus of this line of work, 
however, has been on the optimization of resources (power and 
bandwidth) for practical implementations and on numerical 
computation of resulting capacity bounds (see also fTbl). The 
most relevant to our main question is a recent result by 
Kolte, Ozgiir, and El Gamal ifTTll on a general MIMO relay 
network, which carefully compares the noisy network coding 
lower bound for the general unicast relay network El with 
the cutset bound, which can be readily specialized to the 3- 
node relay channel. In the same vein, another recent study 
by Gerdes, Hellings, Weiland, and Utschick ifTsll establishes 
the optimal input distribution of the partial decode-forward 
lower bound for the MIMO relay channel, the corresponding 
result of which for the single-antenna case is immediate since 
partial decode-forward is the better of decode-forward and 
direct transmission. 

This paper provides more direct and comprehensive answers 
to our main question through an elementary yet careful analy¬ 
sis of the partial decode-forward and compress-forward lower 
bounds for the MIMO relay channel. The main contributions 
are summarized as follows. 

• Eor the complex Gaussian relay channel with ti trans¬ 
mitting antennas at the sender, r 2 receiving and t 2 trans¬ 
mitting antennas at the relay, and receiving antennas 
at the receiver, we show that the partial decode-forward 
achieves within min(fi,r 2 ) bits of the cutset bound 
(Theorem [T]). 

• This gap is somewhat relaxed when noncoherent trans- 
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mission is employed (Proposition | 4 ]i. 

• Unlike the single-antenna counterpart, partial decode¬ 
forward can achieve rates arbitrarily higher than the better 
of decode-forward and direct transmission in MIMO 
relay channels (Proposition [S]). 

• To complement the additive gap result, we show that 
both coherent and noncoherent partial decode-forward 
coding schemes achieve at least half the cutset bound 
(Theorem 121 ). 

• We show that compress-forward achieves min(fi + 

+ r2 bits within the cutset bound (Theorem[ 3 ]). 

• We establish similar results for half-duplex relay channel 
models d, fH, El (Section |Vll). 

In conclusion, the paper establishes simple approximate ex¬ 
pressions of the capacity, which are particularly useful in high 
and low SNR. Beyond these analytical results, we also discuss 
how these expressions can be computed efficiently. 

The rest of the paper is organized as follows. In the next 
section, we formally define the channel model and review the 
cutset upper bound, the partial decode-forward lower bound, 
and the compress-forward lower bound on the capacity. The 
main results on additive and multiplicative gaps for partial 
decode-forward and compress-forward are also stated therein. 
The proofs of these results are given in Sections [III] and |IV] 
Section |V] is devoted to the computational aspects of our 
results, namely, how the capacity bounds can be computed 
efficiently via appropriate convex optimization formulations. 
Using these computational tools, the main results are verified 
by numerical simulations. In Section IVII half-duplex MIMO 
relay channels are discussed. 

Throughout the paper, we use the following notation. The 
superscript (•)'^ denotes the complex conjugate transpose of 
a (complex) matrix; tr(-) denotes the trace of a matrix; 
denotes the nxn identity matrix; denotes a set of nxm 

complex matrices; A> B denotes that A — Bi^, hermitian and 
positive semidefinite; E(-) denotes the expectation with respect 
to the random variables in the argument. 


We assume that the sender and the relay have ti and t2 trans¬ 
mitting antennas, respectively, with average power constraint 
P. As in the standard relay channel model El, the encoder is 
defined by x"(m), the relay encoder is defined by X2i(y2~^), 
i = l,...,n, and the decoder is defined by m{y^). We 
assume that the message M is uniformly distributed over the 
message set. The average probability of error is defined as 
= P{M 7^ M}. A rate R is said to be achievable for 
the relay channel if there exists a sequence of (2"^, n) codes 
such that lim„^oo = 0 . The capacity C of the relay 
channel is the supremum of all achievable rates. 


Z2 



The following upper bound on the capacity is well known. 

Proposition 1 (Cutset bound || 2 ] Th. 4 ]). The capacity C of 
the MIMO relay channel is upper bounded by 

Pcs = sup min{/(Xi,X2; Y3), 

/(Xi;Y 2 ,Y 3 |X 2 )} (2) 

= maxmin{ log 1/^.3 -f \, 

log Kr-2-1-1-3 + G'*iiCi|2G*2 I } (3) 

= imxminj log -f [G31 G32]Pr[G3i G32]^|, 

log |Pi + (G21G21 + G 3 jG 3 i)iCi |2 I} 

( 4 ) 

where the supremum in © is over all joint distributions 
P(xi,X2) such that E(X^Xj) < P, j = 1 , 2 , the maxima 
in and @ are over all (U -|- ^2) X (U + P) matrices 


II. Problem Setup and Main Results 

We model the point-to-point communication system with a 
relay as a MIMO relay channel with sender node 1 , relay 
node 2 , and receiver node 3 ; see Fig. [T] Throughout the 
paper, we assume the complex signal model, but corresponding 
results for the real case can be easily obtained; see the 
conference version ED of the current paper for some results 
on the real model. The relay and the receiver have r2 and r3 
receiving antennas with respective channel outputs 

Y2 =G2iXi+Z2, 

Y3 =G 3 iXi+G 32 X 2 + Z3, 

where G21 G G31 G and G32 G are 

complex channel gain matrices, and Z2 ^ CN( 0 , P2) ^nd 
Z3 ^ CN(0,P3) are independent complex Gaussian noise 
components. For simplicity, we will often use the shorthand 
notation 

G3* = [G31 G32] and G*i = 


G21 

G31 


K = 


Ki 

.-^12 


Ki 2 

K2 


>1 0 


such that tT:{Kj) < P, j = 1 , 2 , and 


Pl|2 = Pi - Pi2P2”^^12- 


( 5 ) 


The equality in ( 0 ]) is justified by the following fact that will 
be used repeatedly throughout the paper. 

Lemma 1 . For 7 G [ 0 , 1 ], r x f matrix G, and t x t matrix 

K to. 


I Ir + jGKG^ \ = \It+jG^GK\ 

>7”in(*’’’)|p + GPGH|. (6) 

We compare the cutset bound with two lower bounds on the 
capacity. The first lower bound is based on the partial decode¬ 
forward coding scheme, in which the relay recovers part of the 
message and forwards it. 
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Proposition 2 (Partial decode-forward bound ^ Th. 7 ]). The 
capacity C of the MIMO relay channel is lower bounded by 

i?PDF = supmin{/(Xi,X2; Y3), 

/(U;Y2|X2) + /(Xi;Y3|X2,U)} ( 7 ) 
= sup min{/(Xi,X2; Y3), 

/(Xi;U,Y3|X2)-/(Xi;U|X2,Y2)} 

(8) 

where the suprema are over all joint distributions F{u, xi, X2) 
such that E(X^Xj) < P, j = 1 , 2. 

Remark 1 . The partial decode-forward lower bound does not 
increase by (coded) time sharing. 

The partial decode-forward lower bound can be relaxed in 
several directions. First, by limiting the input distribution to a 
more practical product form, we obtain the noncoherent partial 
decode-forward lower bound; 


.Rnpdf = supmin{/(Xi,X2; Y3), 

/(U;Y2|X2) + /(Xi;Y 3 |X 2 ,U)} 

( 9 ) 

= supmin{/(Xi,X2; Y3), 

/(Xi;U,Y 3 |X 2 )-/(Xi;U|X 2 ,Y 2 )} 

( 10 ) 

where the suprema are over all product distributions 
E(u,X i)E'(x 2) such that E(XyXj) < P, j = 1 , 2 . Second, 
by setting U = Xi, which is equivalent to having the relay 
recover the entire message, we obtain the decode-forward 
lower bound; 

it-DF = supmm{/(Xi,X2;Y3), /(Xi;Y2IX2)} ( 11 ) 

= maxmin{ log lEg + 

log |-^r 2 + G'2liT3|2G2i I } (12) 

where the supremum in (fTTT i is over all distributions F(xi, X2) 
such that E(X^Xj) < P, j = 1, 2, and the maximum in (fT2l i 
is over all {ti +12) x {ti +^2) matrices /T P 0 of the form Q 
such that ti{Kj) < P, j = 1 , 2 . Third, by setting U = 0 and 
X2 = 0, we obtain the direct-transmission lower bound; 

Rdt = sup J(Xi; Y3) 

= maxlog |/r3 + G3 iP:iG'33| ( 13 ) 

-^1 

where the supremum is over all distributions P(xi) such that 
E(X3 Xi) < P and the maximum is over all ti x ti matrices 
Ki P 0. 

Remark 2 . Since decode-forward and direct transmission 
schemes are two special cases of partial decode-forward, we 
have in general 


-RpDF > max(PDF, .Rdt)- ( 14 ) 

Next, we present another important lower bound, in which 
the relay compresses its noisy observation instead of recover¬ 
ing the message. 


Proposition 3 (Compress-forward bound || 2 l Th. 6], IfTn i. 
The capacity C of the MIMO relay channel is lower bounded 
by 

PCF =sup/(Xi;Y 2 ,Y 3 |X 2 ) ( 15 ) 

where the supremum is over all conditional distributions 
P(xi)P(x2)P(y2|y2,X2) such that E(X^Xj) < P, j = 1,2 
and 

/(X2;Y3)>/(Y2;Y2|X2,Y3). 


This lower bound can be expressed equivalently as 
Rcf = sup min{/(Xi , X2 ; Y3) - /( Y2; Y2 1 Xi , X2, Y3), 
/(Xi;Y2,Y3|X2)} (16) 

where the supremum is over all conditional distributions 
P(xi)P(x2)P(y2|y2,X2) such that E(X^Xj) < P, j = 1 , 2 . 

Remark 3 . The compress-forward lower bound before taking 
the supremum in (flST l or (fTbl l is not a convex function of the 
conditional distribution P(xi)P(x2)P(y2|y2,X2) in general 
and can be potentially improved by (coded) time sharing fT 2 \ 
Remark 16 . 4 ]. 

Remark 4 . By setting Y2 = 0 , compress-forward reduces to 
direct transmission and thus Pcf > ^dt- 

We are now ready to state the main results of the paper. 

Theorem 1 . For every G21, G31, G32, and P, 

ApDF := .Res — ^PDF < niin(fi,r 2 ). (17) 

As a supplement to the additive gap result in Theorem [T] 
which is useful in approximating the capacity in high SNR, we 
establish the following multiplicative gap to provide a tighter 
approximation in low SNR. 


Theorem 2 . For every G21, G31, G32, and P, 

-^^<2. ( 18 ) 
RpDF 

In other words, partial decode-forward always achieves at 
least half the capacity. 


The above results can be relaxed by using the noncoherent 
partial decode-forward. 

Proposition 4 . For every G2i,G3i,G32, and P, 


Anpdf := Res — Rnpdf 

< max[min(fi, r2), min(fi +t2, r3)] 


and 


Res 


< 2 . 


( 19 ) 


( 20 ) 


Rnpdf 

For the single-antenna case, the partial decode-forward 
lower bound can be shown HB Sec. II] to be equal to 
the maximum of the decode-forward and direct-transmission 
lower bounds; cf. (HI. For multiple antennas, however, par¬ 
tial decode-forward is in general much richer than decode¬ 
forward and direct transmission. 


Proposition 5 . If ti,t2,r2,r3 > 2 , 

sup [Rpdf — niax(i?DF, Rdt)] =00. 
G21',Gsi,Gs2,P 
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In IfTTl Th. 1], Kolte, Ozgiir, and El Gamal derived a 
capacity lower bound for general MIMO relay networks. 
When specialized to the three-node relay network and further 
tightened, this result yields the following channel-independent 
capacity approximation. 


Proposition 6. For every G21, G31, G32, and P, 
Acf 

:= Res — Rcf 

ti + t2 


< min max 


min(fi -f f 2 , rs) log 1 -f 


min(fi +t 2 ,r 3 ) 


+ r 2 log(l -f l/cr^), 

min(fi, r 2 -I-r 3 ) log(l-f cr^) . (21) 

We tighten this result further as follows. 

Theorem 3 . For every G21, G31, G32, and P, 

Acf < minmax[min(fi + t 2 , r^) + r 2 log(l + 1 /cr^), 

mm(fi, r 2 -f r 3 )log(l-f cr^)] ( 22 ) 

< min(fi -f f 2 , ^ 3 ) -I- r 2 . (23) 


No multiplicative gap is known between the compress- 
forward lower bound and the cutset bound. This follows partly 
from the fact that the distribution that attains the suprema in 
(fTsT l and (fTbl l is rather difficult to characterize. It can be shown, 
however, that when restricted to Gaussian distributions, the 
compress-forward lower bound (even with time sharing) may 
have an unbounded multiplicative gap from the cutset bound. 
As a compromise, we state the following simple consequence 
of Remark |4] and the proof of Theorem |2] 


Proposition 7 . For every G21, G31, G32, and P, 


Res 

max(i?DF, .Rcf) 


< 2 . 


III. Partial Decode-Forward 


(24) 


In this section, we establish the results on partial decode¬ 
forward stated in the previous section (Theorems [U and |2] and 
Propositions |4] and |5]l. 


A. Partial Decode—Forward (Proof of Theorem\J} 

We evaluate the partial decode-forward lower bound in (| 7 ]l 
with (Xi,X2) ~ CN( 0 , A"), where A" ^ 0 is of the form 
in Q, and 

U = G21X1 + Z'2, ( 25 ) 

where Z2 ~ CN( 0 ,is independent of (Xi, X2, Z2, Z3). 
Note that (Xi,X2,U,Y3) has the same distribution as 
(Xi, X2, Y2, Y3). The first term of the minimum in (| 7 ]i is 

/(Xi,X 2;Y3) =log|R3 +G3*RGS,|. ( 26 ) 

For the second term, since 
Cov(Xi|U,X2) 

= Cov(Xi|Y 2 ,X 2 ) 

= Ari|2 — Ari|2G2i(/r2 + G2lAri|2G2i) ^G2 iAG|2 
= Ri|2(Ri + G2iG2iA'i|2 ^ , 


we have 


/(U;Y2|X2)+/(Xi;Y3|X2,U) 

|R3+G3iCov(Xi|U,X2)GH,| 
IR3 +G2 iCov(Xi|U,X2 )GHJ 
+ log \Ir2 + G2 iKi\2G21 \ 

= log |/ti + (G 21 G 21 + G]^iG3i)Ki\2 I 

. |Ri -I- G2iG2iAri|2| 

^ IRi + 2G21G21AG12I 
4^ log |/ti + (G 21 G 21 + G3]^G3i)Ari|2 I 


(27) 

mm(G,r 2 ) (28) 


where the last inequality follows by Femma[T] Comparing (l26l) 
and ( l28l l with the cutset bound in (HJi completes the proof of 
Theorem [T] 

We can prove Theorem [T] alternatively using the following 
result that is applicable to a more general class of relay chan¬ 
nels and follows by setting p{u\xi,X 2 ) = 2 ^ 2 ) 

in the second form of the partial decode-forward lower bound 
in ®. 

Proposition 8. For a discrete memoryless relay channel 

P{y2,ys\xi,X2) = p{y2\xi,X2)p{y3\xi,X2), 


ApDF := Res — RpDF 

< max I{Xi\U\X 2 ,Y 2 ) 

p{xx,X 2 ) 

where p{u\xi,X2) = PY^\x^,x■iiu\xl,X2)■ 

'How, applying Proposition [ 8 ] to the MIMO relay channel, 
we have 


/(Xi;U|X 2 ,Y 2 ) 

= log|R3 + G21 Cov(Xi IY2, X2)G»i I 
= log|Ri + G2iG2lAri|2(/ti + G2iG2lAri|2) ^ 
_i |Ri + 2G21G21 Ari| 2 | 

|.^il + ^ 21 ^ 2 iA'i| 2 | 

< min(ti, r2). 


B. Noncoherent Partial Decode-Forward (Proof of the First 
Statement of Proposition 0 

We use the following fact. 

Lemma 2 . Let K ^ 0 be of the form in Then, for every 
G31 and G32, we have 


G3iRiGH,+G32R2G^2 AG32R['2G3^+G3iRi2GH2. (29) 


Proof of Lemma \ 2 } Consider 


[G 31 


G32] 


Ki 

— 


-Ki2 

K 2 


[G 31 


G32]^ ^ 0 . 


To prove Proposition 01 let AT ^ 0 be of the form in Q. 
Let Xi ^ CN(0,Ari) and X 2 CN( 0 ,Ar 2 ) be independent 
of each other, and define U as in (l25l l. Then, by Lemmas [T] 
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and| 2 ] the first term of the minimum in the noncoherent partial 
decode-forward lower bound in (|9|l is 


Comparing ( l35T l and ( l36l l with the cutset bound in dUi estab¬ 
lishes that 


/(Xi,X 2 ;Y 3 ) 

= log \Ir3 + G^iKiG^i + G32K2G'^2 I 
> log 


( 30 ) 


+ G 32 i^ 2 G ^2 


+ G32K^2G'ii + G3lifl2GH2) (31) 

> log \Ir 3 + Gs^KG^^ I - min(fi -f t 2 , rs). (32) 

Following steps similar to the coherent case in Section IIII-AI 
we have 

Cov(Xi|U,X 2 ) = Cov(Xi|U) 

= K,{It3 + G^,G2 iKi)-^ 

and 


/(U;Y2|X2) + /(Xi;Y 3 |X 2 ,U) 
J^. 3 +G 3 iCov(Xi|U)G 3 M 
^ +G2 iCov(Xi|U)GH,| 

+ log \ Ir2 + G21K1G21 I 
= log \It 3 + (G^,G2i + G 3 ^G 3 i)iFi| 
. \Ir2 + G21 jFiG2i| 

^ l-^ti + 2 G 23 G 2 iiFi| 

> log + (G^3G21 + G3^G3l)iFl| 


min(fi,r 2 ). (33) 


Comparing (|3^ and (l33l l with the cutset bound in (|4|i com¬ 
pletes the proof. 


C. Multiplicative Gap (Proofs of Theorem\ 2 \and the Second 
Statement of Proposition 0 

By setting U = 0 or Xi in ® and specializing (fTSli to 
independent (Xi,X 2 ), it can be readily checked that i?NPDF 
and max(i?DF, ^dt) are simultaneously lower bounded by 

max I min (log 1/^3 + GsiiFiGgi -f G^ 2 K 2 G^ 2 \^ 

log \ Ir2 + G2liFlG2i I), 

maxlog \Ir3 + GaiATiGgi 11 

Ki ) 


= max mm 
K1.K2 


{ log 17,3 + GaiTFiGg^ + G32772G3"21, 
max(log 17^2 -f G2i7s:iG2i |, 


log 17,.; 


G 3 i 77 ,G 


3 ll)}- 


We further lower bound each term in (l34l i. By OTT i. 

log |7r3 + GsiKiG^i + G‘i2K2G^2 1 


> log 


+ i(G 3 i 77 iGH, + G32772G3H2 
+ G327f5^2G3i + G; 


31AG2G32) 


> ^ log 17,3 + (G3l7fiG3^ + G327f2G3"2 


(34) 


+ G 32 ifr 2 G^i + G 3 l 7 fi 2 G 3 H 2 ) I ■ ( 35 ) 


Similarly, 

max{log|7t3 -f G2 iG2iA:i |, log | 7 t^ +G31G31AGI} 

> i(log |7t3 +G2^G2l7fl| +log \h, +G3^G3iAG|) 

> ilog \lt 3 + (GH,G2i + G3^G3i)AG|. ( 36 ) 


ApDF > max(ANPDF, Adf, Adt) 

> min{ANPDF, inax(ADF, Adt)} > ^-^CS- (37) 

D. Decode—Forward and Direct Transmission (Proof of 
Proposition 0 

Consider the MIMO relay channel with G 31 = diag(( 7 ,1), 
G 21 = diag(l,( 7 ), G 32 = diag(g,g), p > 1 , which is 
equivalent to a product of two mismatched single-antenna 
relay channels, one with the direct channel stronger than the 
sender-to-relay channel and the other in the opposite direction. 
Set Ki \2 = Ki = K 2 = {P/2)l2 in ( |26] ) and ( |28] ), we have 

Apdf > min | log (l -f g '^ P ) (^1 -f (1 -f 5^)^) , 

log(l-f (1-7 5 ^)^) - 2 | 

= log(l + (l + 52)|^'_2. (38) 

In comparison, 

Adf = Adt = max log(l-f Ai)(l-f sf^P 2 ) 
Pl-|-P2<P 

<log(l + A)(l + 52p). (39) 

Therefore, we have 

Apdf - max(ADF, Adt) > log 1 , L/, 2 L “ ^ 

(1 -f A)(l -f 5 ^A) 

which tends to infinity as p —> 00 . Based on this example, 
more examples of larger dimensions can be constructed. 


IV. Compress-Forward 

We prove Theorem [3] Let 77 P 0 be of the form in (|5]l. 
Let Xi ^ CN(0,77i) and X 2 ~ CN( 0 , 772 ) be independent 
of each other, and 

Y 2 = Y 2 + Z 2 (40) 

where Z 2 ~ CN( 0 ,CT^ 7 r 2 ) is independent of Xi, X 2 , Z 2 , and 
Z 3 . Then, 


7 (Y 2 ;Y 2 |Xi,X 2 ,Y 3 ) 

= fi(Y2|Xi,X2,Y3) - fi(Y2|Xi,X2, Y2,Y3) 

= r 2 log(l + l/a2) (41) 


7 (Xi;Y 2 ,Y 3 |X 2 ) 

= log 


(1 + 0 

0 

+ G*i 

TfiGDi 


(1 + ( 7 ^)Ir 2 0 

0 /j-g _ 



= log 


+ (t^^2iG21 + G 3 ^G 3 l)i 7 l 
> log 17*3 + (G^iG 2 i + G 3 ^G 3 l) 77 l | 

- min(fi,r 2 + ra) log(l -f cr^). 


(42) 

(43) 
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The first statement of Theorem [3] is now established by 
substituting ( |3^ . (HTt . and ( l43l l in the compress-forward lower 
bound in (fThl l and comparing it with the cutset bound in 
Setting cr^ = 1 in (|22]) yields the second statement in (l23T l. 


V. Computation of the Capacity Bounds 


A. Formulations of Optimization Problems 


1) Cutset Bound: Computing the cutset upper bound in 
Q can be formulated as the following convex optimization 
problem ifTSll : 

maximize i?cs 


over 
subject to 


where 


Res > 0,iT ^ 0,iTi|2 ^ 0 

i?cs <log|/.3 + G3*i^GH,| 

Res < log|-^r'2+r3 + G*iKi\2G^i\ 
iv{A'lKAi) < P,tr(A^iT.42) < P 
K - AiiTi|2^? ^ 0 


Ih 

0*2 Xtl 


and A 2 


Otl Xi2 
, ^*2 


(44) 


The optimization problem in (l44li can be solved by standard 
convex optimization techniques or packages, e.g., Il23l . 

2) Partial Decode-Forward Lower Bound: Since direct 

computation of © or (H) is intractable, we instead consider 
three lower bounds on i?pDF, namely, i?DF,^DT, and the 
special case of i?pDF evaluated by (l25l l. and take the maximum 
of the three. Note that all three lower bounds can be viewed 
as the partial decode-forward lower bound evaluated by (l25l l 
with a more general choice of Z 2 ^ CN(0, where 

= 00 , 0,1, respectively. Considering more values of 
can further improve the bound at the cost of complexity. 

As for the cutset bound, both i?DF and i?DT can be 
computed efficiently as a convex optimization problem. The 
third bound, characterized by (l26l l and (l27l i. is nonconvex. 
Thus, we evaluate the bound with the optimal solution to 
the convex optimization problem defined by (l26T l and (l28l l. 
A similar approach can be taken for computation of /Jnpdf- 

3) Compress—Forward Lower Bound: We consider two 
convex lower bounds on Rcf, namely, the special case of 
Rcf evaluated by (l40l i with cr^ = 1, namely. 


max min 
Ki,K2 


I log 17^3 + G^iKiG'^i + G327T2G32 1 ~ T2, 
log|7t3 + (iG2^G2l + GH,G3 i)77i|} 


and i?DT (which corresponds to = cxa). As in the case 
of partial decode-forward, considering more values of cr^ can 
further improve the bound at the cost of complexity. 


B. Numerical Results 

We consider the additive and multiplicative gaps on 2000 
2x2 MIMO relay channels with random channel gains 
independently distributed according to CN(0,1). The gaps 
are evaluated by relaxed bounds discussed in the previous 
subsection. The maximum and average of the additive gaps 
are shown in Fig. |2] and similar multiplicative gaps are shown 
in Fig. [3 The simulation results are consistent with the theo¬ 
retical predictions in Theorems [T] |2l and [2 and Proposition 0] 



Fig. 2. The additive gaps between the cutset bound and the partial decode- 
forward and compress-forward lower bounds for randomly generated 2x2 
MIMO relay channels. 



Fig. 3. The multiplicative gaps between the cutset bound and the partial 
decode-forward and compress-forward lower bounds for randomly generated 
2x2 MIMO relay channels. 


VI. Half-Duplex MIMO Relay Channels 

Half-duplex relay channel models are often investigated to 
study wireless communication systems in which relays cannot 
send and receive in the same time slot or frequency band. 
There are two different types of half-duplex models. One 
is the sender frequency-division (SFD) MIMO relay channel 
(Fig. |4(a)| ), in which the channel from the sender to the relay, 
X" Y 2 , is orthogonal to the multiple access channel from 
the sender and the relay to the receiver, (X']^,X 2 ) —?> Y 3 . 
The other is the receiver frequency-division (RFD) MIMO 
relay channel (Fig. |4(b)| i, in which the channel X 2 Yg is 
orthogonal to the broadcast channel Xi (Y 2 , Yg). Both can 
be viewed as special cases of the general (full-duplex) MIMO 
relay channel model. For example, the SFD model follows by 
setting [G31 0 ] G and [O G21] G C’'"xW+*i) 

in ©. Consequently, our main results in Section II continue 
to hold with ti = t'l -b t'l and ra = r-g + Tg for the SFD and 
RFD cases, respectively. 

In the following, we present tighter results that exploit the 
half-duplex channel structure. The proofs are similar to the 
full-duplex case in basic analysis techniques and relegated to 
the Appendix. 
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Z2 



Z2 z^' 



A. Sender Frequency-Division MIMO Relay Channels 

It has been shown by El Gamal and Zahedi ll 24 l that the 
relay channel capacity is achieved by partial decode-forward 
when the sender has orthogonal components. We specialize 
this result to the multiple-antenna case. 

Proposition 9. The capacity of the SFD MIMO relay channel 
is 


C = Res = .RpDF 

sup min{/(X;,X2;Y3), 

F(x'^,X2)F(x'/) 

/(X";Y2)+/(X;;Y3|X2)} 


= max mm 
K 


log 




G 


3 * 


K[ 


K',2 

K. 


iK',2) 

\og\R,+G2iK'^G^,\ 

+ log|/,3+G3iX;|2G3M 


G 3 ^ 


( 45 ) 


where the sup remum is over all F (pcf JC2) F {x.”) such that 

< P and E(X5X2) < P, the 
maximum is over all {t'l + t'l -f ^2) X {t'l + t'i + ^2) matrices 


K[ 

0 

K[2 



0 

A" 

0 

t 0 

( 46 ) 


0 

K2 




such that iT{K[ -f K'() < P, tr(X2) < P, and = K[ — 


We can establish the following (in addition to obvious 
corollaries from the full-duplex results). 

Proposition 10. For every G21, G31, G32, and P, 

Anpdf = G — i?NPDF < min(f3 - 1 -^ 2 , ^3). ( 47 ) 


Proposition 11. For every G21, G31, G32, and P, 

Acf = G - Rcf < min(f'^ + ^2, ^3) -f r-2. ( 48 ) 


B. Receiver Frequency-Division MIMO Relay Channels 
The capacity in this case is not known in general. 

Proposition 12. The capacity C of the RFD MIMO relay 
channel is upper bounded by 

Rcs= sup min{/(Xi;Y')+/(X2;Y"), 

F(xi)F(x2) 

/(Xi;Y2,Y')} 


= max 
K 


ixmin| log IE3 -f G31AGG311 

1 ^ 


-I- max log I/1.3 -I- G32K2G 

K2 


^ I 

32 1 5 


log|/i3 + (G^lG2l + G3^G3l)iTl|} 


( 49 ) 

where the supremum is over all F{xi)F{x2) such that 
E(X^Xj) < P, j = 1 , 2 , and the maxima are over all 
ATi, K2 A 0 such that tr(Aj) < P, j = 1 , 2 . 

As in the cutset bound, coherent transmission is irrelevant. 

Proposition 13. For every G2i,G3i,G32, and P, 

RpDF = Anpdf 


and consequently 

ApDF = Anpdf < min(fi,r2). 

We can further establish the following. 

Proposition 14. For every G2i,G3i,G32, and P, 

Acf < max[min(f 1 , r2 + r'^),r2]. ( 50 ) 


Appendix 

Proof of Proposition 17 ^ Set K'12 = 0 ^nd = K[ 
in (I45I 1. The result follows by similar arguments to (l 3 ^ . ■ 

Proof of Proposition [ 77 } Let K be the form of (l 46 T l. Let 
X'l - CN( 0 ,A:(), X" - CN{ 0 ,K”), and X2 ~ CN(0,A:2) 
be independent and Y2 = Y2 + Z2, where Z2 ~ CN(0 ,/t- 2) 
is independent of (X'^^, X", X2, Z2, Z3). Then, 

/(Xi, X2; Y3) - /(Y2; Y2 1 Xi, X2, Y3) 

= log I A3 + G3iK[G^, + G 32 K 2 G ^2 I - ^2 

and 

/(Xi;Y 2 ,Y 3 |X 2 ) 

= log I A3 + G31 A^G3^ I + log IA2 + (1/2)G2 i A^'G2^ | 

> log I As + G31 a;G 3^ I + log IA2 + G2iK';G»i I 
- min(f", r 2 ). 

The gap due to the first term, min(f']^ -f AjFa) + ^2, follows 
by similar arguments to the proof of Proposition [TOl which 
dominates the gap due to the second term. ■ 

Proof of Proposition IT?} The compress-forward lower 
bound in ([Thi l simplifies to 

Rcf = supmin{/(Xi; Y') + /(X2; Yg)-/(Y2; Y2IX1), 
/(Xi;Y2,Y')} ( 51 ) 

where the supremum is over all conditional distributions 
A(xi)A(x2)A(y2|y2) such that E(X^Xj) < P, j = 1 , 2 . 















Let Xi ^ CN( 0 ,iLi) and X2 ~ CN(0,A"2) be independent 
of each other. Let Y2 = Y2 + Z2, where Z2 ~ CN( 0 , is 
independent of (Xi, X2, Z2, Z3, Z3). Then, 

/(Xi ; Y') + /(X2; Y") - I{Yr,% I Xi) 

= log \Ir 3 + GsiKiG'^i I + log 1/^3 + G' 32 Ar 2 G 32 I ~ ^2 


[22] A. El Gamal and Y.-H. Kim, Network Information Theory. Cambridge: 
Cambridge University Press, 2011. 

[23] M. Grant and S. Boyd, “CVX: Matlab software for disciplined convex 
programming, version 2.0 beta,” http://cvxr.com/cvx Sep. 2013. 

[24] A. El Gamal and S. Zahedi, “Capacity of a class of relay channels with 
orthogonal components,” IEEE Trans. Inf. Theory, vol. 51, no. 5, pp. 
1815-1817, 2005. 


and 


/(Xi;Y2,Y') 



The rest of the proof follows similar steps to that of Proposi¬ 
tion [TT] ■ 
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